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In the past decade we have witnessed remarkable developments in the gauge-gravity duality, 
which suggested a new approach to superstring theory and quantum space-time. In this context 
it is important to study supersymmetric large-A^ gauge theories in the strongly coupled regime. 
I will summarize the results and insights obtained so far by non-lattice simulations. A simple 
example of the gauge-gravity duality is the one between Id U(A^) gauge theory with 16 super- 
charges and the so-called black 0-brane solution in type IIA supergravity. In order for this duality 
to be vaUd, one has to take the 't Hooft large-A' limit and to take the strong coupling limit on the 
gauge theory side. The gauge theory can be regularized by fixing the gauge completely thanks to 
one dimension, and by introducing a Fourier mode cutoff. One can then use the standard RHMC 
algorithm to simulate the system. The energy calculated as a function of the temperature was 
compared with the results obtained from the gravity side based on the black hole thermodynam- 
ics. This confirmed the gauge-gravity duality with high accuracy and provided the microscopic 
origin of the black hole thermodynamics. From the calculation of the Wilson loop, one obtains 
the Schwarzschild radius of the dual geometry. One can actually use the present Id model with 
supersymmetric mass deformation to study — 4 super Yang-Mills theory on RxS^ based on a 
novel large-A^ reduction, which generalizes the original idea of Eguchi and Kawai. A test of this 
approach has been provided by Monte Carlo simulation at weak coupling. It is remarkable that 
we can now simulate the 4d superconformal field theory, which appears in the most typical case 
of the gauge-gravity duality known as the AdS/CFT correspondence. In particular, no fine-tuning 
is required unlike previous proposals based on the lattice regularization. 
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1. Introduction 

I hesitate a bit to talk about non-lattice simulation at this lattice conference, but this is indeed 
crucial for the purpose of simulating supersymmetric theories. What I am going to discuss is the 
so-called gauge-gravity duality, which is a conjecture from superstring theory [|l|, (See ref. 
[^] for a comprehensive review.) The statement itself is simple, and we don't even have to refer 
to superstring theory. Let us consider \]{N) supersymmetric Yang-Mills theory (SYM) with 16 
supercharges (or 32 supercharges, in a special case). We take the so-called 't Hooft limit, which 
amounts to sending N to infinity with fixed X = ^ym^- ^^^"^ consider the strongly coupled 
regime, namely the large-A regime. Then the statement is that the SYM is "dual" to a classical 
solution in lOd supergravity. 

The argument for this conjecture is actually very intuitive and easy to understand. In su- 
perstring theories, there exists a soliton-like object termed D brane. "D" stands for the Dirichlet 
boundary condition imposed at the boundary of the worldsheet of a string. D brane can extend in 
p-\-\ dimensions, and it is characterized as a hypersurface on which strings can end on. Let us con- 
sider an open string attached to the D brane propagating along it. In figure [T] on the left, we describe 
such a process diagrammatically. If one slices the diagram in the orthogonal direction, one notices 
that the same process can be viewed as emission of a closed string. This is an example of the well- 
known notion of open-string/closed-string duality. Note here that an open string includes a gauge 
particle as a massless mode, and similarly a closed string includes a graviton. Let us consider N D 
branes lying on top of each other. Then, in the low energy limit, one obtains {p+ 1) -dimensional 
U(A'^) SYM as an effective theory which describes the massless degrees of freedom of open strings 
attached to the D branes. On the other hand, one obtains a curved lOd space-time in the bulk since 
the D brane sources gravitons. In order for the supergravity to be valid as a low energy and classical 
description of superstring theory in the bulk, one has to take the so-called 't Hooft large-A^ Umit 
with fixed A = g\yiN, and then to take the large-A limit. This is so, since the string loop corrections 
are suppressed by \/N, whereas the a' corrections, which are due to strings having finite extent, 
are suppressed by some powers of 

Why is this duaUty interesting? First of all, it is a realization of an old idea by 't Hooft [Q], 
which states that the large-A^ gauge theory is equivalent to some classical string theory, although 
in those days people may not have anticipated that the string theory actually lives in a curved 
space-time. It is interesting that the curved space-time emerges from a gauge theory in a flat 
space. This aspect of the duality is often referred to as the emergent space-time. In the gauge- 
gravity duality, one typically obtains the anti-de Sitter space. If one considers the gauge theory at 
finite temperature, one obtains a black-hole-like geometry ||5|, ^]. As applications, one can study 
strongly coupled gauge theories, which are relevant to hadron and condensed matter physics, from 
a curved space-time. One can also use the duality in the opposite direction, and try to explain 
the microscopic origin of the black hole thermodynamics in terms of gauge theory. An ultimate 
goal of the gauge-gravity duality is to construct a non-perturbative and background independent 
formulation of superstring theory by using gauge-theory degrees of freedom. 

Since the gauge-gravity duality is a strong-weak duality, it is important to study gauge theories 
in the strongly coupled regime. Monte Carlo simulation can be a powerful tool for such purposes. 
However, the problem is that the gauge theories we are interested in have supersymmetry, which 
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Figure 1: On the left, an open string attached to the D brane is propagating along it. On the right, the same 
process has been viewed as emission of a closed string from the D brane. 



is broken by the lattice. This can be seen immediately if one recalls the supersymmetry algebra 
{Q^Q} °^ P^i^ where the generators for translation appear on the right hand side. Since the transla- 
tional symmetry is broken by the lattice regularization, one necessarily breaks supersymmetry. 

Recently there are considerable developments in "lattice supersymmetry", which can be cat- 
egorized into two classes. One is the construction of lattice actions with various symmetries. For 
instance, one can preserve one supercharge by using the so-called topological twist. (See Catter- 
all's contribution of this volume). The other one, which we discuss here, is non-lattice simulations 
1^] of supersymmetric gauge theories in 1 dimension with 16 supercharges ^ ^, Notably, one 
can extend this approach to 3d and 4d gauge theories JTl] ] by using the idea of reduction 



[12]. In the 4d case, the gauge theory becomes superconformal and the number of supersymmetries 
enhances from 16 to 32. This superconformal theory is interesting on its own right, but it is also 
studied intensively in the context of the AdS/CFT correspondence, which is a typical case of the 
gauge-gravity duality The non-lattice simulation of the 4d superconformal theory requires no 
fine-tuning, unlike the previous proposals based on the lattice regularization [pj|]. 

This article is organized as follows. In section |2| I discuss the non-lattice simulation of Id 
SYM with 16 supercharges. In particular, I explain how black hole thermodynamics appear from 
Id SYM, and how the Schwarzschild radius appears from the Wilson loop. In section ^ I review 
the large-A^ reduction, which enables us to extend these works to higher dimensions. In particular, 
I discuss how one can study ^ = 4 SYM on /? x 5^ in the 't Hooft limit, and present some prelim- 
inary results for the Wilson loop and the two-point correlation functions. In section ^ I conclude 
with a summary. 



2. Non-lattice simulation of Id SYM with 16 supercharges 



The Id SYM with 16 supercharges has the following actions for the bosonic part and the 
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fermionic part, respectively. 

5b = ^ ^'^Jftr|i(DX,(0)'-^[^/(0,^7(0]'} (2-1) 

5f = ^ j\t tv\^^^>aD^>a - \'i'a{yi)ap[Xi,'i'p\^ (2-2) 

It is a Id \i{N) gauge theory, and the covariant derivative is denoted as D = 5, — / [A(?), • ]. 
Xj{t) {] = 1, • • • ,9) and (a = 1, • • • , 16) are N xN Hermitian matrices, and the theory has 

SO(9) symmetry. When we are interested in finite temperature, we impose periodic boundary con- 
ditions onXj{t) and anti-periodic boundary conditions on ^a{t). Then the temperature is given by 
r = j8^\ where f5 is the extent in the Euclidean time (t) direction. The 't Hooft coupling constant 
is defined by A = ^^A'^, which has the dimension of mass cubed. The physics of the system is 
determined only by the dimensionless coupling constant Agff = jj. Therefore one can take A = 1 
without loss of generality. With this convention, the low T regime corresponds to the strongly cou- 
pled regime, which is expected to have the dual gravity description [^], whereas the high T regime 
is essentially weakly coupled, and the high temperature expansion (HTE) is applicable [|T4[]. 

In non-lattice simulation [^, we introduce an upper bound on the Fourier mode as Xi{t) = 
l^^^_^X,-,„e"*^'", where ft) = and similarly for the fermions. This idea does not work usually 
because it breaks gauge invariance. (Recall that the Fourier mode is not a gauge invariant concept.) 
However, in Id, one can fix the gauge non-perturbatively in the following way. We first take the 
static diagonal gauge A{t) = ^diag(ai , • • • , a^), in which the gauge field is constant in time and 
diagonal. By following the usual Faddeev-Popov procedure one obtains 



Sfp = - £ 21n 



a<b 



. Ua-ab 
sm 



(2.3) 



as a term to be added to the action. The above gauge choice does not fix the gauge symmetry 
completely, and there is a residual symmetry given by 

OCa '-^ 0Ca + 2KVa , X"^ i-)- X"^_y^^y^^ , i-)- ^'^^„^v^^v,, ' (2-4) 

which represents a topologically nontrivial gauge transformation corresponding to the gauge func- 
tion g{t) = diag(e"*'^'', • • • ,e'®^"'). This residual gauge symmetry can be fixed by imposing —7i < 
oCa < 71. One can then introduce the Fourier mode cutoff A. Since there is no UV divergence in this 
Id model, one can take the A — oo limit naively, and one obtains the original gauge theory with 16 
supercharges. 

The system with finite A can be simulated efficiently by using the standard RHMC algorithm 
[15]. In particular, the Fourier acceleration [|l^] can be implemented without extra cost since we are 



dealing with the Fourier modes directly as the fundamental degrees of freedom. This is crucial in 
reducing the critical slowing down at large A. (The same theory is also studied using the standard 



lattice approach [17|. However, from the results obtained so far, the non-lattice simulations seem 
to be far more efficient in obtaining the continuum limit.) 

Let us first discuss the phase structure that appears when one changes the temperature. As 
is well known, the Polyakov line serves as an order parameter for the spontaneous breaking of 
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the center symmetry. Figure ^ (Left) shows the results [j8|]. At high temperature the data agree 



nicely with the HTE [|14j] including the next-leading order. As the temperature decreases below 
r ~ 3, the data start to deviate, and at low temperature below T ~ 0.9, the data can be fitted to the 
characteristic behavior of the "deconfined phase" 

a 



(l^l)=exp 



(2.5) 



In the temperature regime investigated, we find no phase transition. This is in sharp contrast to the 
bosonic model [18,1^, |^], which undergoes a phase transition to the "confined phase" at T ~ 0.9. 
The absence of the phase transition is consistent with analyses on the gravity side [0, 21 1. 
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Figure 2: (Left) The Polyakov line is plotted against T . The dashed line represents the result of HTE up 
to the next leading order for = 12 [[l4|]. The dotted line represents a fit to eq. (2^) with a = 0.15 and 
b = 0.072. (Right) The energy (normalized by A^^) is plotted against T. The dashed line represents the result 
obtained by HTE up to the next leading order for = 12 [p^. The solid line represents the asymptotic 
power-law behavior at small T predicted by the gauge-gravity duality. The upper left panel zooms up the 
region, where the power-law behavior sets in. 

Let us turn to a quantitative prediction from the gauge-gravity duality. Given the dual geome- 
try, one can use Hawking's theory of the black hole thermodynamics to obtain various thermody- 



namic relations such as [ ]22[ ] 
E 



1 



A 1/3 



A 1/3 



14/5 



_9_ 
14 



,13152^^ 



14 



1/5 



7.41 ••• 



(2.6) 



The gauge-gravity duality predicts that this should be reproduced by Id SYM in the large-A'^ limit 
at low T The importance of this prediction is that, if it is true, it explains the microscopic origin 
of the black hole thermodynamics, meaning that the Id SYM provides the quantum description of 
the states inside the black hole. 

In figure § (Right) we plot the internal energy |^, which is defined hy E = {^^) in terms of 



the free energy At T > 3 the data agree with the HTE [14]. As one goes to lower temperature. 



the data points approach the solid line, which represents the result (2.6) obtained from the lOd 



black hole. (See refs. [23] for earlier studies based on the Gaussian approximation.) 

The plots in figure ^ were actually presented two years ago at LATTICE 2007 in Regensburg 



[24]. A common criticism in those days was that it was not clear whether the gauge theory results 
continue to follow the line predicted from gravity at lower T . In fact, simulations at lower T 
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are difficult, since one has to increase A proportionally to \/T , and at the same time one has to 
increase to avoid the run-away behavior due to finite N Instead of lowering T, we were able 
to determine the power of the subleading term as [|l( 



1 



A 1/3 



A 1/3 



14/5 



c 



A 1/3 



23/5 



(2.7) 



from gravity. This was derived by considering higher derivative corrections in the supergravity 
action due to the effects of strings having finite extent {a' corrections). The coefficient C of 
the subleading term is calculable in principle, but it requires the full information of the higher 
derivative corrections, which are yet to be determined. By using ( |2.7| ), however, we can already 
make a nontrivial test of the gauge-gravity duality |[Io|]. In figure ^ (Left) we plot the discrepancy 
7.4iri'*/^ —E/N'^ against T in the log-log scale, which reveals that the power of the subleading 
term is indeed consistent with the predicted value 23/5 = 4.6. In figure || (Right) we find that the 
data at r < 0.7 can be nicely fitted to the form ( ^ ) with C = 5.58. Note also that the A = 6 data 
seem to suffer from some finite A effects at low T . From this point of view, we consider that the 
A = 4 data points at low T in figure |^ (Right), which seem to be on the curve of the leading order 
result from gravity, also suffer from finite A effects. Now we know that actually the subleading 



term in (2.7) should be taken into account for precise agreement. 
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Figure 3: (Left) The deviation of the internal energy -^E from the leading term 7.41 is plotted against 
the temperature in the log-log scale for X = \. The solid line represents a fit to a straight line with the 
slope 4.6 predicted from the a' corrections on the gravity side. (Right) The internal energy -^E is plotted 
against T for A = L The solid line represents the leading asymptotic behavior at small T predicted by the 



gauge-gravity duality. The dashed line represents a fit to the behavior (2.7) including the subleading term 
with C = 5.58. 



As another prediction from the gauge-gravity duality, let us consider the Wilson loop, which 
winds around the temporal direction once, like the Polyakov line. However, unlike the usual 
Polyakov line, we consider the one involving the adjoint scalar as 



If dt{A{t) + iniXi{t)] 





(2.8) 



where «, is a unit vector in 9d, which can be chosen arbitrarily due to the S0(9) invariance. This 
object can be calculated on the gravity side by considering the minimal surface spanning the loop 
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in the dual geometry [E 



^. For the present model, the result is given by 



InW 



K 



2na' 

where /?sch is the Schwarzschild radius of the dual black hole geometry and 

2/5 



K - 



1 f 16\/T57r'^/2 V 



271 



7 



7 



li 



(2.9) 



(2.10) 



In figure ^ we plot the log of the Wilson loop [^] against T anticipating (^^. Indeed, at 
low temperature (to the right on the figure), we find that the data points can be fitted nicely to 



a straight line with a slope 1.89 in precise agreement with ( |2.10[ ). The solid line corresponds to 
(log|W|) = 1.89r^^/^ —4.58, where the existence of the constant term can be understood as a' 
corrections. This result demonstrates that one can extract the information of the dual geometry such 



as the Schwarzschild radius from the gauge invariant observable (2.8). It also confirms directly 



the fuzz-ball picture 1 27] of a black hole proposed to solve the information paradox. 



5 

D) 




Figure 4: The plot of (log|W|) for X = \ against T^^l^. The cutoff A is chosen as follows: A = 12 for 
= 4; A = 0.6/r for = 6, 8; A = 4 for = 14; A = 6 for = 17. The dashed line represents the results 
of the HTE up to the next-leading order for A^ =14, which are obtained by applying the method in Ref. [|l4|]. 
The solid line and the dotted line represent fits for N = 6 and N = 4 respectively, to straight lines with the 
slope 1.89 predicted from the gravity side at the leading order. 

One can also predict various correlation functions from gravity. This was done ten years ago 



by Sekino and Yoneya [ 28 1 extending the Gubser-Klebanov-Polykov-Witten prescription M\ to the 



present case. For instance, let us consider an operator 

ffe = Str{Xi,Xi,---Xi,) , (2.11) 

where the symbol "S" implies that all the indices are symmetrized. The two-point correlation 
function of this operator is predicted as 

1 4^-9 
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at A <^\t\ <^X V3/y/^io/2i jj^ figure ^ we plot the two-point correlation function for £ = 4,5, 



which agrees precisely with the predicted power law behavior. See Ref. [ ]29| ] for more details as 
well as results for other operators. 
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Figure 5: The two-point correlation function {ffi{t)ffi{0)) is plotted for £ — 4 (circles) and £ ^5 (triangles) 
in the log-log scale. Simulations were carried out at = 3, A = 12, T = 0.2. The straight lines are fits to 



the predicted power-law behavior ( 2.12 ) 



3. Extension to higher dimensions based on the large-A^ reduction 

In this section we discuss how one can extend the works in the previous section to higher 
dimensions. Respecting supersymmetry becomes more non- trivial in higher dimensions, but here 
again we stick to a non-lattice regularization. For that purpose we use the idea of the \aige-N 
reduction, which we review briefly. Let us consider U(A'^) gauge theory on a D-dimensional torus, 
and consider the Wilson loop defined by 

W[C] = (^^exp(^i J A^{x{o))x''{o)do^'^ , (3.1) 

where the loop C is specified by the embedding function C = {x^ (o)}. The corresponding laige-N 
reduced model can be obtained by simply reducing the torus to a point. This implies that we drop 
the .x-dependence of the field A^(x) and obtain A^. The Wilson loop in the reduced model can be 
defined by 

w[C] = ^^exp (^ij A^x^{o)do^^ . (3.2) 



lim w[Cl = lim IV[C1 . (3.3) 



Then the statement is that 



The original idea was formulated on the lattice by Eguchi and Kawai p(3|]. However, it was soon 
pointed out by Bhanot, Heller and Neuberger pT] ] that there was a problem due to the spontaneous 
breaking of the center symmetry, which invalidates the proof of the statement. Several years ago 



Narayanan and Neuberger [ 32 ] proposed to avoid the spontaneous breaking of the center symmetry 



8 



Non-lattice simulation of supersymmetric gauge theories.. 



Jun Nishimura 



by not reducing the torus to a point completely, but keeping the volume finite in physical units. 
More recently, Kovtun, Unsal and Yaffe pointed out [ ^3| ] that the original proposal with the one- 
site model actually works by adding an adjoint fermion if its mass is sufficiently small, which 



was supported by numerical simulation []34{]. See Bringoltz and Sharpe's contributions as well as 



Hietanen's one on this volume. (See also ref. [ |35| ] for a proposal for any non-abelian gauge theory.) 
Here we use the idea of the reduction in order to study = 4 SYM on R X as 



proposed by Ref. [|12|]. It actually differs from the original large-A'^ reduction in that one deals with 
a curved space rather than a torus, which is a fiat space. The theory obtained after reducing the 
to a point is given by the Id SYM, which is nothing but the one discussed in the previous section, 
plus some mass deformation, which preserves 16 supersymmetries of the undeformed theory. The 
additional terms are given by 



/ 



dt tr 



1 3 J 9 3 



(3.4) 



where is the deformation parameter, which is related to the radius of the before reduction as 

2 
At 



/?S3 = — 



(3.5) 



This mass deformed theory possesses many classical vacua given by X, = /iL,-, where L,- is an 
arbitrary (not necessarily irreducible) representation matrix of the SU(2) algebra \Li^Lj\ = iEijkL^- 
These vacua preserve 16 supersymmetries, and they are all degenerate. 

In order to retrieve the original 4d JV = 4 SYM, one has to pick up a particular vacuum 



(") 



^,=At 



.(n+\) 



- (n+v-1) 



(3.6) 



where lI'"^ represents the m-dimensional irreducible representation of the SU(2) algebra. Note, in 
particular, that there is an identity 



(=1 



(m)\2 



Din 



(3.7) 



which implies that each of L^™^ in (3^) represents a fuzzy sphere with the radius ^\/m?- — I. In 



this construction one regards the as an fiber on S^, where the is represented by a fuzzy 
sphere, and the 5' fibration is represented by having many of them with different radii. A more 



detailed argument for the reduction is given in ref. [|12|]. The statement is that in the /: — 5- oo, ?i — >■ oo 
and V — oo limits, one obtains the ^ = 4 U(oo) SYM on 7? x S^, where the radius of is given by 



( |3.5| ) and the 't Hooft coupling constant is given by 



(2w + v-l)/2 



(3.8) 
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Note that one does not introduce the lattice structure anywhere in the formulation, which is impor- 
tant for preserving supersymmetry. 

A check of this novel Iwge-N reduction has been provided [ jll] ] in the weak coupling limit 
by studying the deconfinement transition at finite temperature. Figure ^ shows that the results 
obtained from the reduced model with the background reproduce the known result [36] for 
^ = 4 U(oo) SYM onRxS^ in the k, n, V — )■ oo Umit. 
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Figure 6: The normalized free energy of the reduced model around the background (3.6) with n ~ ^^-i- in 



the k ^ °o and v ^ °° limits is plotted against the dimensionless parameter x = exp(— ju/2r) representing 
temperature near the critical point Xc = 0.072. The error bars represent the fitting error associated with the 
extrapolation. The solid line represents the result [pq] for the = 4 U(oo) SYM on 7? x S^. 



In order to test the approach at strong coupling, let us consider the circular Wilson loop in 
^ = 4 U(oo) SYM on R^. The expectation value of the Wilson loop is calculated exactly for 



arbitrary coupling constant, and the result is given as [ 37 1 



(^circular) = \j 1 h (V2AsYm) (3.9) 

g-y/2AsYM 

atAsYM>l (3.10) 



(f)V2(2AsYM)^/4 

in terms of the modified Bessel function of the first kind. At strong coupling it agrees with the 
result obtained from the dual geometry [^]. 

Since the ^ = 4 SYM is conformally invariant, the theory on 7?^ is equivalent to the theory 
oaRxS^ through conformal mapping. The circular Wilson loop on is mapped to a great circle 
on at a point on R. (The size of the circular Wilson loop corresponds to the position of the point 
on R, and the dilatation invariance on corresponds to the translational invariance on R.) This 
Wilson loop can be represented in the large-A'^ reduced model in a simple way as 



^circular ]\[^ 



Ak 

exp( /— {X3(0 + /X4(0} 



(3.11) 



where t can be any value due to translational symmetry, and hence one can take an average over it 
to increase statistics. 
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In figure ^ we present our preliminary results for the circular Wilson loop [38]. (Here after, 
we present results obtained by imposing periodic boundary conditions on fermions, since we are 



interested in zero temperature.) We also plot the all order result ([3.9|). Although the matrix size is 
obviously too small, our Monte Carlo results look promising. Note, in particular, that we already 
start to observe a bent from the weak coupling behavior towards the strong coupling behavior. 



A 
V 




Figure 7: The log of the circular Wilson loop normalized by -^/Asym is plotted against V^sym- We have 
performed the A — > oo extrapolation (linear in 1/A) using A = 6,8, 10. The extent in the time direction is 
fixed to j3 =5. The background is chosen to be « = 1, V = 2, A: = 2 for a/Asym < 2, whereas for V^sym > 2, 
we performed an extrapolation Vq k — oa using the data for A; = 2,3 assuming that the finite-A effects are 



0(1 /A ). The solid line represents the all order result (3.9). The dashed line represents the behavior (3.10) 



at strong coupling, whereas the dotted line represents the leading perturbative behavior ln(W) ~ |Asym- 

Next let us consider chiral primary operators such as trZ'', where Z = -^{X\^iX^\ The 
two-point function can be calculated in the weak coupling limit of = 4 U(oo) SYM on 1^ as 



(trZ^(zi)trZt^(z2))«4 



\Z\-Z2 



2J 



Cj 



A. 



SYM 



V 



47r2 ) 



(3.12) 



It is known in ^ = 4 SYM that the supersymmetry non-renormalization theorem holds for the 



two-point functions. Hence, the result (3.12) actually holds for arbitrary coupling constant. As in 
the case of the circular Wilson loop, one can make a conformal mapping to /? x 5^, and obtain for 
7 = 2, for instance. 



^^3 / dQ!^ (trZ2(?,n3)trZt2(0,n'3))«,53 



C2e 



-fit 



(3.13) 



where we have integrated over the since the operator in the reduced model corresponds 

to ^dQ.^tvZ^{t,Q.^) in the SYM on/? x S^. 

In figure |8| we plot our preliminary results for (trZ2(f)trZ^2(0))red against the dimensionless 
time jXt [p9|]. Surprisingly our results for two different values of Asym turn out to be very close to 
each other. This suggests that the non-renormalization theorem actually holds for each background. 
In fact we can obtain results in the weak coupUng limit of the reduced model for the background 
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(3.6), and show that the correlation function approaches the one (3.13) for the ^ = 4 U(oo) SYM 
in the v — )• oo limit ||3^]. Therefore, if we are able to confirm the non-renormalization theorem 
for each background, it immediately implies that we can reproduce the results for the = 4 U(oo) 
SYM from the reduced model at arbitrary 't Hooft coupling constant. 



1 



0.1 



0.01 



0.001 



1e-04 




0.5 



Figure 8: The two-point function of the chiral primary operator is calculated from the reduced model for 



a fixed background (3.6) with n = 3, V = 2 and k — 2. The UV cutoff is chosen to be A = 10. The circles 
represent results for jj. = 3, j5 = 4, which corresponds to Asym = 0.24. The triangles represent results for 
= 1, /3 = 10, which corresponds to Asym = 6.4. Surprising, the results lie more or less on top of each 
other. The solid line represents the analytic result in the weak coupling limit of the reduced model for the 
same background in the A oo and j3 — > 0° limits, which shows reasonable agreement with the Monte Carlo 
data. The dotted line represents the analytic result in the weak coupling limit of the ^ — 4 U(°°) SYM, 
which is expected to be reproduced from the reduced model in the v — > 0° limit. 

What we have presented so far should be considered as a check of our method. More interest- 
ing quantities are those which are not obtained in the strongly coupled gauge theory, and yet there 
exist interesting predictions from gravity. Calculating such quantities by our method will clearly 
provide a new test of the AdS/CFT correspondence. 

For instance, we can consider the rectangular Wilson loop in R^, which behaves as 

{W{TxR))=exp(^^^ (3.14) 

at r » /? due to conformal symmetry. This is in striking contrast to the area law in pure Yang-Mills 
theory. In particular, the AdS/CFT correspondence predicts 

47rV2AsYM 



r4(l/4) 



(3.15) 



at strong coupling [25] 



It would be also interesting to study higher point functions of the chiral primary operators. (See 



ref. [40] for calculation of extremal 3-point functions by simulating a "truncated theory" composed 
of six commuting bosonic matrices.) In particular, AdS/CFT predicts that non-extremal 4-point 
functions violate the non-renormalization theorem. It is interesting to check whether this is indeed 
the case, and to obtain explicit results, which can be compared with the prediction from gravity. 
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4. Summary 

I hope I have convinced the readers that non-lattice simulations are indeed useful for studying 
supersymmetric large-A'^ gauge theories in the strongly coupled regime. In the first part, I dis- 
cussed the Id SYM with 16 supercharges, which reproduced black hole thermodynamics and the 
Schwarzschild radius of the dual geometry. These results revealed direct connections to gravity 
as predicted by the gauge-gravity duaUty. In particular, the gauge theory results provided a clear 
understanding of the microscopic origin of the black hole thermodynamics. In the second part, 1 
discussed how one can extend these works to higher dimensions by using the novel large-A^ re- 
duction. I have presented some prehminary results for ^ = 4 U(oo) SYM on RxS^. This theory 
is superconformal, and it actually has 32 supersymmetries. Our formulation preserves 16 super- 
symmetries in the A — > oo Umit, and the remaining half of the supersymmetries are expected to be 
restored without fine-tuning by increasing the matrix size. 

In a way, what we have seen is the beginning of a whole new field of research analogous to the 
situation of the lattice gauge theory in early 80s. Now with the aid of supersymmetry and large A^, 
we have just started to explore superstring theory and quantum space-time from first principles. 
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